Suppression of the D'yakonov-PereP spin relaxation mechanism for all spin 
components in [111] zincblende quantum wells 
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We apply the D'yakonov-Perel' (DP) formaUsm to [lll]-grown zincblende quantum wells (QWs) 
to compute the spin lifetimes of electrons in the two-dimensional electron gas. We account for both 
bulk and structural inversion asymmetry (Rashba) effects. We see that, under certain conditions, 
the spin splitting vanishes to first order in k, which effectively suppresses the DP spin relaxation 
mechanism for all spin components. We predict extended spin lifetimes as a result, giving rise to the 
possibility of enhanced spin storage. We also study [110]-grown QWs, where the efi'ect of structural 
inversion asymmetry is to augment the spin relaxation rate of the component perpendicular to the 
well. We derive analytical expressions for the spin lifetime tensor and its proper axes, and see that 
they are dependent on the relative magnitude of the BIA- and SIA-induced splittings. 

PACS numbers: 85.75.-d,73.21.Fg 
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I. INTRODUCTION 

If the current pace of electronic device miniaturization 
is to continue, it is reasonable to think that the good 
use of the quantum properties of the electron will play 
a role in making this possible. Traditionally it has been 
the wave character of the electron that has been put to 
this use, resulting in devices such as the resonant tunnel 
diode^ and the single electron transistor—. 

Another quantum property of the electron that only re- 
cently has received attention for its potential for informa- 
tion storage and processing is its spin. The study of spin 
dynamics in two-dimensional electron gases (2DEGs) is 
crucial to the understanding of the operation of spin 
electronic (spintronic) devices such as the Datta-Das 
transistor—. A newly proposed family of devices^ is 
based on the special properties of the spin hfetime tensor 
due to the interplay between bulk inversion asymmetry^ 
(BIA) and structural inversion asymmetrj^ (SIA) in a 
[001] quantum well (QW), as pointed out by Averkiev 
and Golub^, and Kiselev and KimS. The effects of SIA 
on the spin dynamics should always be kept in mind, as 
it can be unintentionally present in any heterostructure 
due to uneven doping profilesiS^ surface effects, differ- 
ent intcrdiffusion at the boundaries, etc. Previous stud- 
ies have found analytical expressions for the D'yakonov- 
Perel' (DP) spin lifetimesii in quantum wells (QWs) tak- 
ing into account BIA effects onljsi^, or considering both 



BIA and SIA for [001] QWs only^. There have also 
been numerical-only studies of the spin lifetimes in [110] 
QW£i3. Experimental studies in [110] structures have 
shown evidence of extended spin lifetimes^^ and lifetime 
modulation by the application of a gate bias^^. 

In this work we investigate how the interplay between 
SIA and BIA affects the spin lifetimes for electrons in 
[111] and [110] QWs. We start by computing the ef- 
fective spin Hamiltonians in a two-band model for the 
[111] and the [110] cases. We then proceed to compute 
the ensemble lifetime of the three spin components as a 
function of the relative magnitude of BIA and SIA con- 
tributions following the procedure from Refs. IT3 andlTR 
We finally discuss the results and the device implications 
of our findings. 



II. EFFECTIVE SPIN HAMILTONIANS 

To find the effective two-band Hamiltonians for the 
zincblende QWs, we start from the 0{k^) spin part of 
the Hamiltonian for bulk zincblendes^^ 



BIA 



7 [(^xkx [kl - kf) + c.p] , 



(1) 



where ai are the Pauli matrices, ki are the electron 
wavevector components and c.p. stands for the cyclic per- 
mutation needed to obtain the remaining terms of the 
Hamiltonian. 
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TABLE I: Correspondence between growth orientation- 
dependent X, y, z labels and crystallographic orientations. 



We first do a change of basis to express Hbia in natu- 
ral coordinates for the [111]- and [110]-grown structures. 
Then, following the procedure in Refs. and 0, we 
quantize k along the growth direction and, keeping only 
terms linear in k|| — second order terms in ky vanish be- 
cause of time reversal requirements for the expectation 
value of — , we arrive at the following expressions for 
the BIA Hamiltonian of [111] QWs 



BIA [1111 



and [110] QWi 



M. 



BIA [1101 



V3 



{ky<7x kxCTy) , 



(2) 



(3) 



where the labels x, y, z depend on the orientation of the 
structure (see Table QJ. 

Upon inspection of the Hamiltonians in Eqs. (0)-® 
we see that BIA causes a k-dependent effective magnetic 
field pointing in-plane for [111] structures, while it points 
along the growth direction for [110] structures. Note that 
Hbia [hi] is formally identical to the Rashba Hamilto- 



(4) 



where aR is the Rashba coefficient, whose value depends 
on the particulars of the structural asymmetry present 
in the sample. We shall now see that this has important 
consequences in the values of the spin lifetimes. 



III. [lll]-GROWN STRUCTURES 

We define aeiA = 27(fcf)/-\/3. At this point we drop 
the [110] and [111] subindices where it is clear from the 
context to which structure we refer. The combination of 
Eqs. (|2J) and Q yields the first order Hamiltonian for 
[111] structures 



HiA,i = (obia + Or) {kyUx ~ kxOy) 

= aiA {ky(T.j: - kxUy) , 



(5) 



where we have introduced a parameter ajA = olbia + ctR 
describing the combined effects of BIA and SIA in the 
heterostructure. 



Since Eq. ^ is formally identical to the Rashba Hamil- 
tonian, all the results that have been derived for the case 
of SIA onljiSiii will also hold even when BIA is accounted 
for, just by making the substitution ctR 



Q!ia: 



1 



2a2 kHi ' 



(6) 



where the tilde indicates a magnitude that is evaluated at 
a given energy and f i is the effective time for field reversal 
due to the harmonic / = 1 of the scattering cross section. 



and in genera. 



17.21 ~-i 



[E) = §a{4>,E){l - cosl<t>)d<t> 



The spin directions will be perpendicular to the wavevec- 
tor and in-plan o^^i^^ . 

Doing the thermal average of the corresponding scat- 
tering rates^^ and assuming fi oc we obtain the spin 
lifetimes for a nondegenerate 2DEG population. 



= = 2t2 = 



1 



e 'p 



(^ + 1)! 



n> 1, 



(7) 



(8) 



with Tp = {ly + 1)! fi(kg) being the transport time, the 
symbol ( )g meaning a (non-degenerate) thermal av- 
erage, fcg = 2m*kBT/h^ is a thermal wavevector, T is 
the electron temperature, ks is the Boltzmann factor, 
and m* is the electron effective mass. For the degener- 
ate case the thermal averaging is trivial and one obtains 
Ti — Ti{kp), with kp being the Fermi wavevector. We 
see that, as usual in the DP mechanism, the spin lifetime 
is inversely proportional to the momentum lifetime. 

A most interesting configuration for [lll]-grown sam- 
ples occurs when q;bia = — Q;r,. Then, aiA = and the 
conduction bands become spin degenerate to first order 
in k. The most significant consequence of this configura- 
tion is that the spin lifetimes would be extended for any 
spin direction, as opposed to spins along [110] for (100) 
structures and Q!bia — ctR^ , or spins perpendicular to the 
plane well for (110) structures and or = Oi^. Control of 
q;r can be achieved by the application of a gate bia322i24 
or by sample design with compositional asymmetry, pro- 
viding a nonzero an at zero bias. Thus, properly biased 
(111) QWs could act as spin reservoirs, or form the basis 
of a resonant spin lifetime transistor as has been recently 
proposed for [100}M and [110] structures''^. 

We have performed 8-band k • p calculations that show 
the validity of our treatment. Figure ^ shows the con- 
duction band spin splitting as a function of the angle 
of k|| with respect to k^ for zero applied perpendicular 
bias and -0.38 V, where we find that aeiA ~ — ctR- R is 
clearly seen that forV — OY {i.e. or = 0) the splitting is 
fairly independent of the angle, an indication of the pre- 
dominance of the linear regime. On the other hand, for 
V — —0.38 V we see a much smaller splitting (note the 
logarithmic scale) with an oscillatory character indicative 
of higher order in k terms dominating the Hamiltonian 
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FIG. 1: a) CB splitting as a function of angle at = 
0.005 A*^ for an 8ML AlSb/GaSb/AlSb quantum well. 
Shown are results for V and -0.38 V of applied bias (note 
the logarithmic scale), b) Value of oia as a function the per- 
pendicular bias. 



[see Eq. @ below]. The inset shows the calculated de- 
pendence of aiA (through or) on the voltage applied, 
illustrating the control of ajA with a gate voltage. We 
believe that the high values of aeiA are due to [111] being 
the direction of maximum asymmetry in the zincblende 
structure. 

We can find the limitation that higher order terms im- 
pose upon these resonant spin lifetimes by following the 
procedure of Averkiev tt al^^. For (111) structures, the 
effective CB Hamiltonian will be 



V2 (3fc2 _ fc2) fc^^^ 



(9) 



If we add Eq. ® to Eq. we obtain the following 
results for the spin lifetimes 



6^2 



1 



fc2fi I2a2^ _ 4V37aiAfc2 + (1 + 2^3/^1 )7^fc^ 
3ft2 1 



fc2fi (^^2 _ 2V3aiA) 



(10) 
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FIG. 2: Spin lifetimes for the three spin components for [001]-, 
[110]- and [lll]-grown QWs as a function of the ratio of the 
SIA and the BIA parameters. Note the resonance for all com- 
ponents for in the case of a [111] structure, as opposed to 
a single component resonance for [001] and [110] structures. 
The labels r^^n^c, for the [110] orientation refer to the spin 
relaxation tensor proper axes as defined in Eq. 1131 . 



Since the scattering rate is proportional to {Hia)'^, the 
terms in Eq. (|1U|) are not correct in general because 
terms arising from the combination of ffiA,i with fifth 
order contributions to Hia are missing. However, we 
have kept the terms here because they are correct in 
the special case where aiA = 0, giving the lowest order 



4 



contribution to the spin scattering rate. 

Figure 121 shows the calculated lifetimes in the present 
treatment for the two structures under study, and previ- 
ous [001] results^ are provided for convenience. Equa- 
tion (|10|l is plotted in Fig. |2]a) for typical values kp = 
0.01 A-\ Tp = 1 ps, 7 = 186 eV-A^^, ubia = H x IQ-^" 
eV-ciEr^ as a function of the ratio aR/aeiA- The three 
spin lifetime components show the predicted resonant 
spin lifetime when aR = — qibia- Eq. (|10|l by itself limits 
the lifetimes at the resonances to finite but very large val- 
ues. In practice, other mechanismsSLSSiSS will effectively 
limit the value of the resonant spin lifetime. Therefore, 
[lll]-grown hctcrostructures provide DP suppression on 
par with [110]-grown structures, with the added advan- 
tage the suppression is for all spin components, as op- 
posed to one component only. 



IV. [110]-GROWN STRUCTURES 

We also apply the present formalism to derive analyt- 
ical expressions for the spin lifetimes in [110] QWs un- 
der the combined action of BIA and SIA. Here we define 
obia = and the effective Hamiltonian due to 

inversion asymmetry is 



which induces a spin splitting Aja.i equal to 



A 



IA,1 



= 2k 



"BIA 



COS^ 



(11) 



(12) 



where k and are defined by k^ = k cos and ky = 
ksiiKf). 

Figure 13 shows the spin splitting and directions for the 
lower conduction band of a [110]-grown QW for different 
values of r = aR/aeiA and k fixed to 0.01 A~^. In plot 
a), the distance from the line to the center indicates the 
amount of splitting, and the angle corresponds to (j>. For 
r = we see that spin splitting is highly anisotropic, com- 
ing from the BIA contribution of Eq. lfTT|l . As the SIA 
effects, which are inherently isotropic, grow, the splitting 
becomes more and more isotropic. The spin directions 
for r = [plot b)] have very interesting behavior. The 
effective magnetic field always points along ±z, indepen- 
dently of kH. This will result in absence of precession for 
spins along z, and the DP mechanism for spin relaxation 
will be suppressed for that component of spin^S. 

The proper axes of the spin lifetime tensor are depen- 
dent on the value of or and q;bia- Their orientation i 
and their corresponding lifetimes fi are found by diago- 
nalizing the spin scattering rate tensor: 
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R + Q^BIA, 



aBiAU + aRZ 
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(13) 
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FIG. 3: Spin splitting and directions for the lower conduction 
subband of a [110]-grown QW. Plot a) shows the value of the 
spin splitting from Eq. Ill2t for states lying on a circle in the 
kx, ~ ky plane with k = 0.01 for various values of the 

ratio QfR/aBiA. Plots b) and c) show the spin eigenstates for 
Eq. im for OR/aBiA = and aa/aBiA = 0.35, respectively. 
Note that when an, — all spins are aligned with the z axis. 
In that case, no spin direction is specified for the [010] and 
[010] directions because the states are spin degenerate. The 
value assumed for aeiA is 10~^ eV-cm. 



Equation (|13() is plotted in Fig. |2|b) for the same val- 
ues as in Fig. [21a). The and components show 
regular behavior, while the t,^ component shows the res- 
onant spin lifetime when or = 0i2i2S that has been jus- 
tified above. If an is due to some external bias, we see 
from Eq. (|13|) that the spin lifetime along direction ( will 
decrease proportionally to E^'^, where E is the applied 
electric field perpendicular to the QW plane. In practice, 
this resonance will have a finite peak height due to the 
other mechanismsSLSSiSr starting to kick in. In principle, 
higher order in k terms might also limit the spin lifetime 
at the resonance. However, we see that the k^ terms that 
should be added to Eq. ifTTl) are given by 



H'L = 7 (-kl/2 + kl) k. 



(14) 



'BIA 



The direction of the effective magnetic field is still z ir- 
respective of the value of k. As a consequence, the reso- 
nance is not destroyed by the inclusion of the k^ terms, as 
noted by Hall et alm^. As for the case of [111] structures, 
the possible effect of higher than k^ terms in the Hamil- 
tonian would be shadowed by the onset of the other spin 
relaxation mechanisms. 
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V. SUMMARY 

In conclusion, we have calculated and given analyti- 
cal expressions for the electron spin lifetime for [111] and 
[110] QWs including both BIA and SIA (Rashba) effects. 
We find that the D'yakonov-Perel' (DP) spin relaxation 
mechanism can be effectively suppressed for all spin com- 
ponents in [111] QWs at a resonance condition through 
appropriate sample design or the application of a suit- 
able gate bias. The z spin component of [110] QWs can 
be made to show the same level of suppression of the 
DP mechanism, but the other two components will al- 
ways have the non-resonant dynamics. This effect in [111] 
QWs can lead to the relaxation of design constraints in 



devices such as the spin-LED or the resonant spin life- 
time transistor. Also, these structures may act as a spin 
storage stage. 
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